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ABSTRACT 


The  angular  distribution  of  dose  albedo  was  measured  for  "semi- 
infinite"  concrete,  aluminum  and  steel  slabs  irradiated  by  plane- 
parallel  beams  of  Co°®  and  Csl37  gamma  photons .  The  photons  were 
incident  on  the  slabs  at  angles,  measured  from  a  normal  to  the  slab,  of 
arccos  1.00,  0.75  and  0.50.  For  the  necessary  sensitivity,  the  reflected 
dose  rate  was  measured  by  a  digital  dosimetry  system  using  a  plastic 
scintillator  as  the  detector. 

An  empirical  formula  for  angular  dose  albedo  was  derived  from  the 
experimental  data.  Comparisons  are  made  with  the  Monte  Carlo  calcu¬ 
lations  of  D.  J.  Raso  and  the  experimented  results  of  C.  E.  Clifford 
(DRCL). 
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SUMMARY 


In  many  shielding  problems  it  is  important  to  know  the  amount  of 
gamma  radiation  "reflected"  (backscattered)  from  a  surface  relative  to 
that  which  is  incident  upon  the  surface.  Ihis  quantity,  when  measured 
in  terms  of  dose,  is  referred  to  as  dose  albedo.  In  this  esqperiment 
dose  albedos  have  been  measured  for  thick  slabs  of  concrete,  aluminum 
and  steel  at  various  angles  relative  to  the  perpendicular  to  the  slab's 
surface.  Two  gamma-ray  sources,  Co*3^  and  Os' 37,  were  each  used  to 
irradiate  the  slabs  at  three  angles  of  incidence. 

In  order  to  facilitate  the  solution  of  those  shielding  problems 
which  involve  backscattered  gamma  radiation,  a  simple  mathematical 
expression  has  been  derived  from  the  experimental  data.  The  experimental 
data  have  also  been  compared  with  theoretical  calculations  that  vise  the 
Monte  Carlo  technique.  In  general,  the  experimental  results  are  higher 
than  the  theoretical  calculations  and,  in  some  cases,  exceed  the 
theoretical  values  by  more  than  50  percent. 
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I.  INTRODUCTION 


This  investigation  is  one  of  a  series  arising  from  the  basic 
experimental  shielding  program  at  NRDL.  Previous  studies  have  been 
concerned  with  transmitted  gamma  photon  spectral  and  angular  distri¬ 
bution,  ^  with  transmitted  and  reflected  dose, 2  and  with  the  angular 
distribution  of  transmitted  dose. 3 

This  report  presents  experimental  measurements  of  the  angular 
distribution  of  dose  reflection  coefficients  (albedos)  for  "semi¬ 
infinite"  concrete,  aluminum  and  steel  slabs  irradiated  by  plane- 
parallel,  monoenergetic  gamma  photons.  The  sources,  Co“°  and  Cs^-37, 
were  each  collimated  to  irradiate  nearly  all  of  the  surface  of  a  slab. 

A  plastic  scintillator  was  collimated  to  point  at  the  center  of  a  slab 
face  and  to  "view"  an  area  of  a  slab  surface,  which  was  small  with 
respect  to  the  irradiated  surface  area. 

Since  one  of  the  objectives  of  this  investigation  was  to  provide 
an  experimental  verification  of  Raso*s  Monte  Carlo  calculations,^  the 
polar  angles,  0,  and  azimuth  angles,  0,  for  detector  positioning  and 
incidence  angles,  0O,  for  source  positioning  were  chosen  to  correspond 
to  those  used  by  Raso . 

An  exact  source  energy  correspondence  was  not  feasible  for  Raso's 
choice  of  input  energies.  However,  Chilton  and  Huddleston, 5  losing  a 
semi -empirical  formula  and  Raso's  data,  have  calculated  dose  albedos 
for  Co°*  and  Csl37  source  energies.*  These  albedo  values  for  concrete 
have  been  tabulated  in  this  report  for  comparison. 

The  format  of  this  report  was  adopted  as  an  attempt  to  lead  the 
reader  to  the  results  as  quickly  as  possible.  Details  concerning  the 
means  by  which  these  results  were  derived  are  relegated  to  the 
appendices.  Appendix  I  contains  the  notation  and  definitions  used  in 
this  report  as  well  as  a  discussion  of  the  concept  of  gamma-ray  albedo. 
In  Appendix  II  is  found  a  discussion  of  the  conditions  imposed  upon  the 


*  Private  communication. 
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design  of  the  experimental  apparatus  by  the  albedo  quantity  to  be 
determined,  and  the  methods  by  which  they  were  satisfied. 


II.  DOSIMETOY  SYSTEM 


The  inputs  to  the  dosimetry  system  are  the  reflected  gamma-ray 
photons  from  the  slab  face.  The  system  is  composed  of  four  sections: 
detector,  converter,  register,  and  reductor.  The  function  of  the 
detector  section  is  to  convert  the  energy  deposited  in  the  plastic 
scintillator  by  each  interacting  photon  into  an  electrical  pulse  whose 
voltage  amplitude  is  proportional  to  that  energy.  The  function  of  the 
converter  section  is  to  amplify  and  convert  these  voltage  pulses  into 
digital  pulses.  In  the  conversion  process  a  pulse  train  is  generated, 
the  number  of  pulses  in  the  train  being  a  digital  representation  of  the 
deposited  energy.  Each  conversion  also  gives  rise  to  a  pulse  which  is 
called  an  events  pulse.  In  the  register  section  the  digital  and  events 
pulses  are  accumulated  in  the  digital  and  events  registers,  respectively. 
In  the  reductor  section  these  accumulated  quantities  are  combined  to 
determine  detected  dose  rate,  which  is  designated  dD  to  conform  with 
Raso's  notation. 


III.  DATA  PRESENTATION 

A.  Calculation  of  C^(fi) 


The  quantity  to  be  determined  from  the  experimentally  measured 
quantities  is  the  differential  angular  dose  albedo,  Q&(fi).  %(fi)  is 

defined*  as  the  ratio  of  Ja(fi),  the  differential  angular  "dose -current" 
scattered  ("reflected")  out  of  the  slab,  to  (J^)  ,  the  "dose -current" 
incident  on  the  slab.  In  equation  forcn : 


*  See  Discussion  of  Albedo,  Appendix  I. 
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CtAa)  -  Jd^  keV/g.s.sr 
(Jd;o  keV/g.s 


The  quantities  measured  were  dD,  the  detected  dose  rate,  and  (l^)0,  the 
dose  flux  incident  on  the  slab.  It  is  now  necessary  to  establish  the 
relationships  between  the  measured  and  desired  quantities.  (jd)Q  is 
readily  obtained  from  (la)0  by  the  relationship  between  current  and 
flux.  Therefore 


(jd>0  =  (ld>oCOS6o 


where  6  is  the  angle  between  the  source  to  slab  center  line  and  the 
slab  normal  as  shown  in  Fig.  1. 

In  order  to  relate  dD  to  Ja(ft)  the  following  is  noted:  the 
radiation  received  by  the  detector  is  limited  to  that  which  is  within 
a  small  but  finite  solid  angle,  ft6  (effective  detector  solid  angle), 
defined  by  the  detector  collimator;  further,  Ia(ft)*  is  the  angular 
dose  rate  scattered  from  the  slab  and  I^ftjdft  is  the  dose  rate  from 
the  center  of  the  slab  in  direction  ft,  within  the  solid  angle  d 
Therefore,  the  following  relationship  exists  between  dD  and  I^(ft). 

ft 

€ 

dD  =  /  Id(ft)dft. 

o 

Now,  assuming  that  I^(n)  is  constant  in  the  integration  interval, 


dD  =  id(n)n€. 


*  See  Discussion  of  Albedo,  Appendix  I. 
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Finally,  using  the  relationship  between  flux  and  current  quantities, 


dD 


cose 


n 


e 


) 


where  @  is  the  angle  between  the  detector  to  slab  center  line  and  the 
slab  normal  as  shown  in  Fig.  1. 

Therefore,  the  desired  relationship  is 


dD  cose 

n  (i.)  cose  ’ 

€  d'o  o 


A  more  ^"physical"  quantity  used  in  the  analysis  of  the  experimental 
data  is  A^(fi).  As  stated  by  Raso,^  "...  A<j*  is  the  dose  rate  ratio 
[angular  dose  reflection  coefficient]  which  would  be  physically  observed 
by  §n  infinitesimal  directional  detector  embedded  in  the  surface . " 

Afl(fl)  is  defined  as  follows: 


cosQ  ^ 

cosG2-  “d^* 


therefore 


Ad(n) 


dD 

"  fl#  Cid)0‘ 


*  Raso  uses  A^  and  Q&  where  this  report  uses  A^(n)  and  Ck(fl).  Our 
notation  is  intended  to  alert  the  reader  to  the  fact  that  the  values 
are  per  steradian.  Wherever  values  derived  from  Raso  are  quoted,  his 
notation  will  be  used. 
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B.  Tabulated  Data 


Tables  1  through  4  and  7  through  10  contain  values  of  dose 
albedo  and  Tables  5  and  6  contain  values  of  dose -rate  ratio.  Hie  values 
are  displayed  as  entries  in  arrays  in  which  the  column  headings  are  the 
detector  ©  values  and  the  row  headings  are  the  detector  0  values.  Each 
table  contains  three  arrays-,  one  for  each  source  incidence  angle,  eo. 

The  arrays  for  ©0  =  0  (normal  incidence)  have  only  one  0  row  since  in 
this  case  azimuthal  symmetry  prevails.  Hie  experimental  values  have 
been  multiplied  by  ICr  and  the  calculated  values  by  10^  for  presentation 
in  the  tables. 

Tables  1  and  2  contain  oa(fl)  values  for  12  in.  of  concrete 
obtained  from  the  experimental  data.  Tables  3  and  4  contain  o&  values 
predicted  by  Chilton-Huddleston  formula  based  upon  Raso 's  Monte  Carlo 
calculation  for  a  semi -infinite  concrete  medium.  Hie  Chilton-Huddleston 
formula  is 


ck(©8)  1026  +  C' 

°a  =  1  +  cos©  sec© 
o 


where  K(©s)  is  the  KLein-Nishina  differential  energy  scattering  cross 
section  for  the  scattering  angle  ©s.  C  and  C*  are  fitting  parameters. 

In  Tables  1,  3,  5,  7,  9  the  monoenergetic  source  photon  energy 
is  662  keV  ( 08^37 ) ,  in  Tables  2,  6,  8  and  10  the  source  was  Co°°  with 
photon  energies  1170  and  1330  keV.  In  Table  4  the  Chilton-Huddleston 
formula  used  1250  keV  as  an  equivalent  monoenergetic  source  photon 
energy. 


In  Tables  5  and  6  are  displayed  the  %(fi)  values  corresponding 
to  the  o &(C)  values  in  Tables  1  and  2.  In  Tables  7  through  10  experi¬ 
mental  values,  0&(n),  are  presented  for  aluminum  and  steel  slabs  for 
selected  incident  and  detection  angles .  Hie  underlined  values  are  the 
experimental  values.  Those  values  not  underlined  are  calculated  from  an 
empirical  formula  derived  in  Section  IVA. 
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Table  1 


DOSE  ALBEDO,  a,(0) 
k  d 
(xlO4) 


Source: 

Cs13? 

Material : 

Concrete  (p*155 

lb/ft3) 

Energy: 

662  keV 

Thickness : 

12  in. 

cos©  =  1 

O 

* 

,.00 

Jf\e  _ 

8° 

i?° 

22° 

32° 

45° 

60° 

72° 

7-5° 

** 

NT 

NT 

NT 

156 

137 

102 

66 

cos©  =  0 
0 

'.75 

8° 

15° 

22° 

32° 

45° 

60° 

72° 

7-5° 

221 

220 

223 

229 

234 

220 

182 

22.5 

222 

221 

222 

222 

220 

210 

169 

37-5° 

222 

219 

220 

216 

212 

190 

148 

52.5° 

220 

219 

2l4 

210 

201 

172 

127 

67-5° 

217 

213 

210 

201 

189 

154 

no 

32.5° 

215 

214 

208 

200 

181 

142 

NT 

97.5° 

216 

209 

203 

193 

NT 

NT 

NT 

112 -5° 

2l4 

209 

NT 

NT 

NT 

NT 

NT 

127-5° 

211 

NT 

NT 

NT 

NT 

NT 

NT 

i42.5° 

212 

NT 

NT 

NT 

NT 

NT 

NT 

157-5° 

212 

NT 

NT 

NT 

NT 

NT 

NT 

172.5 

210 

NT 

NT 

NT 

NT 

NT 

NT 

cos©  =  0 
0 

.50 

8° 

15° 

22° 

32° 

45° 

60° 

72° 

7.5° 

297 

305 

324 

356 

402 

502 

516 

22.5 

295 

302 

321 

347 

370 

442 

435 

37.5° 

293 

299 

306 

324 

34o 

362 

332 

52.5° 

291 

296 

299 

304 

307 

308 

236 

67-5° 

288 

290 

283 

284 

272 

244 

182 

82.5° 

286 

285 

268 

269 

244 

206 

145 

97.5° 

280 

275 

262 

256 

226 

183 

126 

112.5° 

278 

269 

259 

248 

213 

NT 

NT 

127.5° 

272 

2  68 

262 

240 

NT 

NT 

NT 

142.5° 

270 

263 

259 

234 

NT 

NT 

NT 

157-5° 

265 

263 

255 

228 

NT 

NT 

NT 

172.5 

265 

258 

253 

226 

NT 

NT 

NT 

*  For  normal  incidence,  albedo  is  independent  of  0. 

**  NT  =  data  not  taken  because  of  interference  between  source  and 
detector  or  supports. 
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Table  2 

DOSE  ALBEDO.  aJsl) 
k  d 
(xiO4) 

Source:  Co^°  Material:  Concrete  (p=155  lb/ft^) 

Energy:  1170  and  1330  keV _ Thickness:  12  in.  _ 


* 

cose  =  1.00 

0 

Ae  ... 

8° 

15° 

22°  32° 

45° 

6o° 

72° 

7.5° 

#* 

NT 

NT 

NT  99 

89 

68 

42 

C036  =0.75 

0 

8° 

— 

22°  32° 

45° 

6o° 

0 

CVJ 

7.5° 

151 

154 

157 

164 

162 

I67 

145 

22-5° 

151 

153 

156 

l6l 

155 

159 

128 

37-5° 

149 

154 

153 

155 

153 

142 

111 

52.5° 

146 

150 

148 

l48 

144 

127 

93 

67.5° 

146 

146 

144 

l4o 

132 

110 

79 

82.5° 

144 

142 

l4o 

134 

123 

97 

NT 

97-5° 

143 

139 

137 

127 

NT 

NT 

NT 

112-5° 

14l 

137 

NT 

NT 

NT 

NT 

NT 

127-5° 

l4l 

NT 

NT 

NT 

NT 

NT 

NT 

142.5 

l4o 

NT 

NT 

NT 

NT 

NT 

NT 

157-5° 

l4i 

NT 

NT 

NT 

NT 

NT 

NT 

172.5 

139 

NT 

NT 

NT 

NT 

NT 

NT 

cose  =  0.50 
o 


8° 

15° 

22° 

32° 

45° 

6o° 

72° 

7.5° 

219 

234 

255 

289 

345 

44o 

48o 

22.5° 

215 

234 

248 

280 

314 

389 

393 

37-5° 

214 

228 

236 

255 

291 

308 

297 

52.5° 

211 

219 

226 

234 

249 

24l 

207 

67.5° 

204 

211 

212 

214 

210 

193 

143 

82.5° 

198 

202 

203 

188 

182 

152 

109 

97.5° 

193 

195 

195 

175 

162 

131 

88 

112.5° 

194 

191 

184 

167 

148 

NT 

NT 

127-5° 

191 

185 

179 

164 

NT 

NT 

NT 

1.42.5° 

192 

183 

172 

158 

NT 

NT 

NT 

1-57.5° 

187 

180 

174 

156 

NT 

NT 

NT 

172.5 

I89 

180 

171 

150 

NT 

NT 

NT 

*  For  normal  incidence,  albedo  is  independent  of  0. 

**  NT  =  data  not  taken  because  of  interference  between  source  and 
detector  supports. 
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Table  3 


# 

CALCULATED  DOSE  ALBEDO,  O, 

c  d 

137  (xlO  ) 

Source:  Cs  Material:  Concrete 

Energy :  662  keV  _ Thickness :  oo 


cos© 

0 

=  1.00 

8° 

lf>° 

22° 

320 

45° 

60° 

72° 

O0 

7.5° 

1491 

1476 

1452 

cosG 

0 

1399 

=  0.75 

1289 

1079 

8o4 

349 

8° 

15° 

22° 

32° 

45° 

60° 

72° 

84° 

7.5° 

1788 

1801 

l8ll 

1821 

1829 

1817 

1682 

980 

22.5 

1786 

1797 

l8o4 

1807 

1800 

1754 

1584 

894 

37.5° 

1782 

1789 

1790 

1782 

1749 

1652 

1431 

764 

52.5° 

1777 

1778 

1773 

1750 

1690 

1542 

1276 

64o 

67.5° 

1771 

1766 

1753 

1717 

1632 

1444 

1149 

545 

32.5° 

1764 

1753 

1733 

1686 

1582 

1367 

1056 

1*81 

97.5° 

1758 

1741 

1715 

1658 

1541 

1310 

993 

440 

112.5 

1752 

1730 

1699 

1636 

1510 

1270 

950 

4i4 

127-5° 

1747 

1721 

1 686 

1618 

i486 

1241 

923 

398 

142, .5° 

1742 

1713 

1676 

l6o4 

1469 

1222 

905 

388 

157-5° 

1740 

1709 

1669 

1596 

1459 

1211 

894 

382 

172.5 

1738 

1706 

1666 

cos© 

0 

1591 

=  0.50 

1453 

1205 

889 

379 

8° 

15° 

22° 

32° 

45° 

6o° 

72° 

CO 

7.5° 

2220 

2283 

2360 

2512 

2831 

3413 

3777 

2579 

22.5 

2215 

2271 

2337 

2463 

2711 

3123 

3295 

2145 

37.5° 

2205 

2249 

2297 

2380 

2521 

2699 

2633 

1581 

52.5° 

2192 

2220 

2246 

2282 

2321 

2301 

2067 

1135 

67-5° 

2177 

2188 

2193 

2189 

2150 

2005 

1684 

858 

82.5° 

2162 

2157 

2144 

2108 

2020 

1808 

1455 

704 

97-5° 

2146 

2128 

2101 

2044 

1927 

1684 

1322 

621 

112.5 

2132 

2102 

2065 

1995 

1862 

1606 

1245 

576 

127.5° 

2120 

2082 

2038 

1959 

1818 

1556 

1199 

550 

i42.5° 

2111 

2067 

2018 

1934 

1788 

1525 

1171 

535 

157«5° 

2105 

2057 

2004 

1918 

1770 

1506 

1154 

527 

172.5 

2102 

2051 

1998 

1910 

1761 

1498 

1147 

523 

*  Monte  Carlo  calculations  of  Raso  interpolated  by  Chilton-Huddleston 
semi -empirical  formula  with  C  =  0.0390  C*  =  0.0170. 

**  For  normal  incidence,  albedo  is  independent  of  0. 
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Table  4 


CALCULATED  DOSE  ALBEDO,  a  * 

60  ( xio5 )  11 

Source :  Cow  v  '  Material:  Concrete 

Energy:  1250  keV _  _ Thickness ;  oo _ 

cos©  =  1.00 
o 


8° 

15° 

22° 

32° 

45° 

60° 

72° 

84° 

7.5° 

833 

827 

8l8 

796 

749 

648 

500 

227 

cos© 

0 

=  0.75 

8° 

15° 

22° 

320 

45° 

/-  0 

60 

72° 

00 

'c6 

7.5° 

1047 

1072 

1098 

1138 

1197 

1264 

124o 

782 

22.5° 

1045 

1067 

1090 

1124 

1170 

1209 

1151 

696 

37.5° 

104i 

1059 

1076 

1098 

1122 

1118 

1016 

574 

52.5  ° 

1035 

1047 

1057 

1066 

1065 

1018 

879 

462 

67.5° 

1028 

1033 

1035 

1031 

1008 

927 

765 

377 

82.5° 

1021 

1019 

1013 

997 

956 

852 

679 

320 

97.5° 

1014 

1005 

993 

967 

912 

794 

617 

282 

112 *5° 

1007 

993 

975 

94i 

876 

751 

574 

256 

3.27.5° 

1001 

982 

960 

920 

850 

719 

543 

24o 

142.5° 

996 

973 

948 

905 

830 

697 

523 

229 

157.5° 

993 

968 

94o 

895 

818 

684 

511 

222 

172. 5° 

991 

965 

936 

890 

812 

677 

505 

219 

cos© 

0 

=  0.50 

8° 

15° 

22° 

32° 

45° 

6o° 

72° 

%■ 

CO 

7.5° 

1364 

1434 

1520 

1681 

2001 

2615 

3196 

2538 

22.5 

1359 

1422 

1498 

1637 

1896 

2344 

2680 

1961 

37.5 

1349 

i4oo 

1459 

1560 

1728 

1959 

2022 

1307 

52.5° 

1336 

1372 

l4io 

1469 

1549 

1607 

1501 

862 

«7-5° 

1320 

1339 

1357 

1379 

1392 

1343 

1163 

610 

S2.5° 

1303 

1306 

1306 

1298 

1266 

1162 

956 

474 

97.5° 

1287 

1275 

1261 

1231 

1171 

io4o 

831 

398 

H2.5° 

1272 

1248 

1222 

1179 

1102 

958 

752 

354 

3.27.5° 

1259 

1226 

1192 

1138 

1053 

904 

702 

326 

142.5° 

1249 

1209 

1169 

ino 

1019 

868 

670 

309 

157-5° 

1242 

1198 

1154 

1092 

998 

847 

651 

300 

172.5° 

1238 

1192 

1147 

1083 

988 

836 

642 

295 

*  Monte  Carlo  calculations  of  Raso  interpolated  by  Chilton-Huddleston 
semi -empirical  formula  vith  C  =  0.0630  and  c'  =  0.0093. 

**  For  normal  incidence,  albedo  is  independent  of  0. 
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Table  5 


DOSE-RATE  RATIO,  A,(fi) 
L  d 

(xlO  ) 


Source • 

Cs13? 

Material : 

Concrete  (p=155  lb/ft3 

Energy: 

662  keV 

Thickness : 

12  in. 

cose  = 
0 

1.00* 

8° 

15° 

22° 

32° 

45° 

60° 

72° 

7-5° 

#* 

NT 

NT 

NT 

184 

193 

204 

213 

cose  = 
0 

0.75 

8° 

15° 

22° 

32° 

45° 

6o° 

72° 

7-5° 

168 

171 

181 

202 

249 

330 

44i 

22.5° 

168 

171 

180 

197 

233 

314 

4li 

37.5° 

168 

170 

178 

191 

225 

284 

359 

52.5° 

166 

170 

174 

186 

213 

258 

307 

67.5° 

l64 

167 

170 

178 

200 

231 

267 

62.5° 

163 

166 

168 

177 

193 

213 

NT 

97.5° 

164 

163 

164 

171 

NT 

NT 

NT 

112.5° 

162 

163 

NT 

NT 

NT 

NT 

NT 

127-5° 

160 

NT 

NT 

NT 

NT 

NT 

NT 

1.42.5° 

161 

NT 

NT 

NT 

NT 

NT 

NT 

157-5° 

161 

NT 

NT 

NT 

NT 

NT 

NT 

172.5 

159 

NT 

NT 

NT 

NT 

NT 

NT 

cose  * 
0 

0.50 

8° 

15° 

22° 

32° 

45° 

60° 

72° 

7.5° 

150 

158 

175 

210 

284 

502 

835 

22.5 

149 

156 

173 

204 

262 

442 

703 

37.5° 

148 

155 

165 

191 

24o 

362 

537 

52.5° 

147 

153 

161 

179 

217 

308 

382 

67.5° 

145 

150 

152 

167 

193 

244 

294 

82.5° 

144 

147 

145 

159 

173 

206 

235 

97*5 

142 

142 

l4i 

151 

160 

183 

20^ 

112.5 

l4i 

139 

l4o 

146 

150 

NT 

NT 

1.27.5° 

138 

139 

l4i 

142 

l44 

NT 

NT 

i42.5° 

136 

136 

l4o 

138 

NT 

NT 

NT 

3-57.5° 

134 

136 

138 

134 

NT 

NT 

NT 

172.5 

134 

135 

136 

133 

NT 

NT 

NT 

*  For  normal  incidence,  albedo  is  independent  of  <f>. 

**  NT  =  data  not  taken  because  of  interference  between  source  and 
detector  supports . 
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Table  6 


DOSE-RATE  RATIO,  A  (n) 

(xlO4) 

5o  o 

Source:  Co  Material:  Concrete  (p=155  Ib/ff3) 

Energy:  1170  and-  1330  keV _ Thickness:  12  in. _ 


cose  = 
o 

1.00* 

0\e 

8° 

15° 

22° 

32° 

45° 

6o° 

72° 

7.5° 

** 

NT 

NT 

NT 

cosG  = 
0 

117 

0.75 

126 

136 

137 

8° 

15° 

22° 

32° 

45° 

6o° 

72° 

7.5° 

ll4 

120 

127 

145 

172 

251 

351 

22.5° 

ilk 

119 

126 

143 

165 

238 

312 

37.5° 

113 

119 

124 

137 

162 

213 

271 

52.5° 

110 

ll6 

120 

131 

153 

190 

225 

67.5° 

ill 

113 

ll6 

124 

i4o 

165 

191 

82.5° 

109 

111 

113 

119 

130 

i46 

NT 

97.5° 

109 

108 

111 

113 

NT 

NT 

NT 

U2.5° 

107 

107 

NT 

NT 

NT 

ITT 

NT 

3.27.5° 

107 

NT 

NT 

NT 

NT 

NT 

NT 

^2.5° 

106 

NT 

NT 

NT 

NT 

NT 

NT 

3-57.5° 

106 

NT 

NT 

NT 

NT 

NT 

NT 

172.5° 

105 

NT 

NT 

NT 

NT 

NT 

NT 

cose  =  0.50 

O 


8° 

15° 

22° 

32° 

45° 

6o° 

72° 

7.5° 

111 

121 

137 

170 

244 

44o 

776 

22.5 

109 

121 

134 

165 

222 

389 

637 

37.5° 

108 

118 

127 

150 

20  6 

308 

48i 

52.5° 

106 

ll4 

122 

138 

176 

241 

335 

67-5° 

103 

109 

115 

126 

14 Q 

193 

231 

82.5° 

100 

105 

110 

111 

12 9 

152 

176 

97-5° 

98 

101 

105 

103 

115 

131 

142 

112.5° 

98 

99 

99 

98 

104 

NT 

NT 

3.27.5° 

96 

96 

96 

97 

NT 

NT 

NT 

3-42.5° 

96 

94 

93 

93 

NT 

NT 

NT 

3-57*5° 

94 

93 

94 

92 

NT 

NT 

NT 

172.5 

95 

93 

92 

89 

NT 

NT 

NT 

*  For  normal  incidence,  albedo  is  independent  of  0. 

*#  NT  =  data  not  taken  because  of  interference  between  source  and 
detector  supports. 
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Table  7 

DOSE  ALBEDO,  0,(0) 
k  d 
(xlO4) 

(underlined  values  are  experimental  values) 

137 

Source:  Cs  Material:  Aluminum 

Energy:  662  keV  Thickness:  12  in. 

* 

cos©  =  1.00 
o 


0  \® 

8° 

15° 

22° 

32° 

45° 

6o° 

72° 

7.5° 

185 

181 

175 

162 

137 

101 

65 

cos©  = 
0 

0.75 

8° 

15° 

22° 

32° 

45° 

6o° 

720 

7-5° 

235 

235 

227 

225 

231 

226 

221 

190 

22.5° 

227 

224 

227 

221 

211 

178 

37.5° 

234 

225 

223 

223 

212 

196 

155 

52.5° 

232 

224 

220 

219 

201 

179 

132 

67-5° 

232 

223 

217 

215 

191 

163 

116 

32.5° 

231 

220 

215 

210 

186 

152 

103 

97.5° 

231 

219 

211 

207 

181 

145 

95 

112.5° 

231 

219 

211 

205 

177 

l4i 

91 

127 .5° 

231 

219 

211 

205 

176 

138 

89 

i42.5° 

231 

218 

210 

202 

174 

138 

87 

157-5° 

231 

218 

209 

200 

173 

137 

87 

172.5 

231 

216 

208 

199 

172 

136 

87 

cos©  = 
0 

0.50 

8° 

15° 

22° 

32° 

45° 

6o° 

72° 

7.5° 

297 

308 

302 

346 

413 

485 

502 

22.5° 

295 

305 

299 

332 

380 

422 

m 

37.5° 

293 

299 

284 

312 

335 

333 

305 

52.5° 

293 

292 

273 

287 

293 

254 

211 

67.5° 

289 

286 

260 

265 

252 

198 

150 

82.5° 

287 

280 

248 

246 

223 

160 

115 

97-5° 

285 

276 

243 

236 

206 

137 

95 

H2.5° 

28l 

272 

236 

227 

197 

128 

86 

127-5° 

279 

269 

232 

220 

191 

121 

82 

142.5° 

277 

265 

228 

219 

190 

119 

80 

3-57.5° 

277 

263 

226 

217 

187 

117 

78 

172.5 

277 

26l 

226 

217 

184 

115 

77 

*  For  normal  incidence,  albedo  is  independent  of  0. 
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Table  8 

DOSE  AIBEDO,  0,(0) 
j.  d 
(xlO4) 

(underlined  values  are  experimental  values) 

Source:  Co  Material:  Aluminum 


Energy: 

1170  and  1330  keV 

Biickness : 

12  in. 

COS©  = 
0 

1.00* 

8° 

15° 

22° 

32° 

45° 

60° 

72° 

7-5° 

109 

106 

104 

cos©  = 
0 

97 

0.75 

82 

62 

42 

8° 

15° 

22° 

32° 

45° 

6o° 

72° 

7-5° 

144 

155 

152 

163 

178 

171 

146 

22.5° 

144 

155 

151 

159 

172 

1ST 

154 

37.5° 

143 

153 

150 

155 

154 

145 

108 

52-5° 

l4i 

152 

147 

152 

153 

128 

89 

£7-5° 

l4i 

151 

145 

147 

143 

113 

72 

82.5° 

l4o 

148 

142 

142 

133 

101 

60 

97-5° 

l4o 

147 

138 

139 

133 

95 

52 

H2.5° 

i4o 

147 

138 

137 

129 

91 

47 

!27.5° 

i4o 

147 

138 

137 

128 

87 

45 

U2.5° 

l4o 

146 

137 

135 

126 

87 

43 

157-5° 

l4o 

l46 

136 

132 

124 

86 

43 

172.5° 

l4o 

144 

135 

cos©  = 
0 

131 

0.50 

124 

85 

43 

8° 

15° 

22° 

32° 

45° 

0 

0 

72° 

7*5° 

21 6 

232 

250 

282 

342 

450 

476 

22.5° 

214 

230 

247 

268 

310 

3BE 

355 

37*5° 

212 

224 

232 

248 

254 

299 

272 

52.5° 

212 

216 

221 

222 

222 

220 

179 

£7-5° 

208 

211 

208 

200 

181 

164 

117 

82.5° 

206 

205 

197 

181 

153 

126 

82 

97*5° 

204 

201 

191 

171 

136 

103 

62 

112.5° 

200 

195 

183 

163 

126 

94 

53 

127.5° 

198 

193 

ISO 

156 

120 

87 

49 

^2.5° 

196 

189 

176 

154 

119 

85 

47 

157-5° 

196 

187 

174 

153 

116 

83 

46 

172.5 

196 

185 

174 

153 

113 

81 

44 

*  For  normal  incidence,  albedo  is  independent  of  0. 
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Table  9 

DOSE  ALBEDO,  0,(0) 
k  d 
(xlO4) 

(underlined  values  are  experimental  values) 

Source:  Cs1^  Material:  Steel 

Energy :  662  keV _ Thidmess:  $.9  in. 

cos©  =  1.00* 


0\© 

8° 

15° 

22° 

320 

45° 

6o° 

72° 

7.5° 

135 

131 

128 

110 

28 

81 

50 

cos©  = 
0 

0.75 

8° 

15° 

22° 

320 

45° 

6o° 

72° 

7-5° 

1 66 

170 

164 

162 

185 

189 

m 

22.5° 

T55 

170 

163 

158 

180 

179 

159 

57.5° 

165 

169 

162 

154 

171 

163 

136 

52.5° 

164 

167 

159 

150 

159 

I5S 

114 

67.5° 

164 

166 

157 

l46 

150 

131 

97 

82.5° 

162 

164 

155 

l4i 

144 

119 

85 

97*5° 

162 

162 

151 

138 

l4o 

113 

77 

112.5° 

162 

162 

150 

136 

136 

109 

73 

127-5° 

162 

162 

150 

136 

135 

105 

70 

3.42.5° 

162 

161 

150 

133 

133 

105 

68 

157-5° 

162 

161 

148 

131 

131 

io4 

68 

172.5 

162 

160 

147 

130 

130 

103 

68 

COS©  = 

0 

0.50 

8° 

15° 

22° 

32° 

45° 

6o° 

72° 

7.5° 

224 

240 

261 

299 

345 

471 

46i 

22.5° 

222 

235 

258 

285 

312 

4o8 

375 

57.5° 

220 

232 

243 

265 

2S7 

319 

258 

52.5° 

220 

224 

232 

239 

225 

2W 

164 

67.5° 

216 

218 

219 

217 

183 

184 

103 

82.5° 

212 

213 

208 

198 

156 

146 

68 

97-5° 

210 

209 

202 

188 

139 

123 

48 

112.5° 

208 

203 

194 

180 

129 

114 

39 

127.5° 

206 

201 

191 

173 

123 

107 

35 

142. 5° 

204 

197 

187 

172 

122 

105 

33 

157*5° 

204 

195 

185 

170 

119 

103 

32 

172.5 

204 

193 

185 

170 

116 

101 

30 

*  For  normal  incidence,  albedo  is  independent  of  0. 
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Table  10 


DOSE  ALBEDO.  a,(fi) 
h  d 
(xlO4) 

(underlined  values  sire  experimental  values) 

Source :  Co^  Material:  Steel 


Energy: 

1170  and  1330  keV 

Thickness : 

5.5  in, 

cose  = 
0 

# 

1.00 

8° 

15° 

22° 

32° 

45° 

6o° 

72° 

7*5° 

78 

76 

75 

cose  = 
0 

§1 

0.75 

58 

47 

32 

8° 

15° 

22° 

32° 

45° 

6o° 

72° 

7.5° 

120 

125 

117 

122 

145 

159 

148 

22.5 

120 

125 

116 

119 

139 

149 

136 

37.5° 

119 

124 

115 
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IV.  DATA  ANALYSIS 


A.  Concrete 

Since  a  model  of  albedo  has  been  suggested'*  which  utilizes 
single  Compton  scattering  as  one  of  the  important  processes,  it  seemed 
natural  to  attempt  to  plot  the  experimental  data  as  a  function  0B.  0g 
is  the  angle  appearing  in  the  Compton  formula 


\  =  \  =  —  (l-cose  ). 

o  me  v  s' 


From  Fig.  1  it  may  be  seen  that  ©s  is  related  to  0Q,  9,  and  0  by  the 
following  formula: 


cos©s  m  sin©osin9cos0  -  cos0ocose. 


Various  plots  of  the  data  in  terms  of  (%(n)  gid  A^fi)  versus 
0S  and  cos0s  were  tried.  It  was  found  that,  when  A^(fi)  was  plotted 
versus  0S,  all  of  the  points  for  the  same  source  energy,  ho0,  and 
incidence  angle,  0Q,  fell  near  a  single  smooth  curve. 

Further  investigation  yielded  an  empirical  formula  for  these 

curves : 


Ad(fl)  =*  c  e“me®  +  b. 


(1) 
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Since  the  form  of  this  equation  was  suggested  by  a  similar  dependence 
upon  scattering  angle  of  transmitted  angular  dose*  for  0  =  0  the 
parameters  c,  m  and  b  were  evaluated  for  each  (hoQ, ©0 )  set  by  fitting 
these  points  for  which  0  =  7*5  <leg. 

In  Table  11  the  values  found  for  these  parameters  are  shown. 
The  solid  curve  in  each  of  Figs.  2-5  has  been  calculated  using  formula 
(l)  and  the  appropriate  parametric  values  from  Table  11. 


Table  11 

EMPIRICAL-FIT  PARAMETERS 
TO  EXPERIMENTAL  Ad  VALUES  FOR  CONCRETE 


cose 

c 

m 

b 

o 

Cs13? 

r  i.oo 

.102 

1.660 

.0171 

1  -75 

.330 

2.131 

.0142 

l  .50 

I.260 

3.362 

.0133 

Co60  j 

1 

.202 

1.970 

.0105 

agm 

Ml 

2M2 

.0098 

mm 

.739 

2.892 

.0085 

The  exponential  variation  with  0g  is  most  strikingly 
demonstrated  by  plotting  [a<i( fi) -bl  versus  0S  in  radians  (Figs.  6-9). 
In  each  of  these  figures  the  b  level,  which  has  been  subtracted,  is 
indicated  as  a  dashed  line.  From  our  experiences  with  the  dosimetry 
system  we  attribute  the  scattering  of  points  around  the  e“m0s  line  to 
instrumental  instability  rather  than  to  physical  relationships. 


*  Unpublished  investigation  by  N.  E.  Scofield. 
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Fig.  2 


Aa(fl)  vs.  9S.  The  experimental  data  are  plotted  with  symbols 
corresponding  to  the  9  angles  at  which  they  were  taken. 


0.09 


Pig.  3  Ad(fl)  vs.  ©s .  The  experimental  data  are  plotted  with  symbols 
corresponding  to  the  e  angles  at  which  they  were  taken. 


0 


2< 


NRDL  379-64 


Fig.  6  [a^( vs.  08.  The  experimental  data  are  plotted  with  symbols 
corresponding  to  the  0  angles  at  which  they  were  taken. 
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NRDL  379-64 


Fig.  7  [Afl(n)-b]  vs.  9S.  The  experimental  data  are  plotted  with 

symbols  corresponding  to  the  9  angles  at  which  they  were  taken. 
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NRDL  379-64 


Fig.  8  vs.  6S.  The  experimental  data  are  plotted  with 

symbols  corresponding  to  the  0  angles  at  which  they  were  taken. 
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9%  (RADIANS) 


Fig.  9  vs •  ®g*  The  experimental  data  are  plotted  with  symbols 

corresponding  to  the  6  angles  at  which  they  were  taken. 
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B.  Aluminum  and  Steel 


Measurements  using  concrete  as  the  scattering  medium  were 
completed  and  analyzed  prior  to  any  measurements  using  aluminum  or 
steel.  Thus  the  empirical  relationship  between  the  scattering  angle 
and  dose-rate  ratio  had  already  been  noted  and  the  necessity  for 
measurements  at  more  than  representative  space-angle  (@8 )  positions  was 
obviated.  The  positions  selected  were  those  for  which  underlined  values 
of  ctd(fl)  appear  in  Tables  7  through  10. 

In  the  process  of  deriving  the  c,  m  and  b  parameters  for  these 
aluminum  and  steel  data  It  was  observed  that  the  shape  of  their  curves 
approximated  those  of  any  set  of  concrete  data  points.  This  suggested 
that,  If  suitable  average  values  c  and  m  were  selected,  all  of  the  dose 
albedo  data  points  taken  could  be  fit  by  a  single-parameter  formula: 


c  e 


’s  +  b 


where  now  only  b  is  a  function  of  E  ,  0  and  Z.  This  simplification  la, 
necessarily,  accompanied  by  some  lo8s  1 8  accuracy.  In  Pig.  10  the 
differences,  [Ad(n)-b],  between  the  measured  dose -albedo  values  and  the 
fitting  parameter,  b,  are  plotted  as  points  and  the  function,  e”*^,  is 
drawn  as  a  solid  curve.  The  empirically  determined  values  of  b  for  each 
(£q,90j_Z)  configuration  are  tabulated  In  the  figure.  The  selection  of 
1  for  c  and  «  for  m  was  made  by  an  approximation  process .  Ho  fundamental 
significance  Is  to  be  Inferred  from  their  presence  in  the  formula  since 
they  were  chosen  partly  for  their  mnemonic  value. 

Although  no  extensive  effcjrt  was  made  to  optimize  the  fit, 

90  percent  of  the  data  points,  A^fl),  fall  within  15  percent  of  the 
value,  A&(tt),  predicted  by  the  formula  for  the  same  angle,  Qg . 

Thus  the  values  of  Q^(n)  In  Tables  7  through  10  not  experi¬ 
mentally  measured  have  been  calculated  by  the  formula: 


%(<>) 


r 


cose 

cose 

O 


e-*®s  +  b 

. 


where  the  values  of  b  for  each  combination  of  source,  scattering  medium, 
V  and  e  were  determined  from  the  experimentally  measured  value  for  that 
combination. 
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V.  ERROR  ESTIMATION  AND  DISCUSSION  OF  RESULTB 


A.  Systematic  Error 

The  probable  systematic  error  in  albedo  depends  in  part  on  the 
probable  error  in  the  experimental  (14)0*8.  The  strength  of  the  sources 
vas  determined  by  measuring  the  dose  rate  at  various  distances  from  the 
sources  using  calibrated  Landsverk  and  Baldwin-Farmer  dosimeters .  These 
measurements  had  a  probable  error  of  0.3  percent,  but  it  is  unlikely 
that  the  source  strengths  were  known  on  an  absolute  basis  to  better 
than  2  percent.  An  even  greater  source  of  systematic  error  is  in  Kg, 
the  system  dose  conversion  factor.  The  method  of  finding  Kjj  and  its 
probable  error  is  described  in  Appendix  II,  D.  The  probable  error  of 
1.4  percent  refers  only  to  K(hu, T)^  for  123. 5  keV  to  662  keV  photons. 

The  degree  to  which  the  reductor  formula  predicts  dose  for  photons 
outside  of  this  energy  region  is  not  known.  However,  from  previous 
spectral  measurements^-  for  space  angles  up  to  ©B  =  6o°,  it  1b  known 
that  the  photon  number  flux  density  for  scattered  radiation  falls  off 
rapidly  for  photon  energies  below  100  keV;  further,  all  scattered 
photons  from  Csl37  and  most  scattered  from  Co°°  for  the  experimental  08 
values  have  energies  less  than  662  keV.  Thus,  even  if  KChn^T)  should 
vary  considerably  from  the  observed  trend,  we  estimate  the  overall 
probable  systematic  error  to  be  no  greater  than  10  percent. 


B .  Relative  Error 


Many  dose -rate  measurements  were  repeated  during  the  course  of 
the  experiment.  The  repetitions  were  made  after  various  combinations 
of  source  and  detector  repositioning  and  gain  adjustment  after  equip¬ 
ment  breakdown  or  time  lapse.  The  measurements  were  found  to  be 
reproducible  to  within  1  percent. However,  if  one  considers  the 
scatter  of  data  points  in  the  A^(0)  versus  0O  presentations  of  Figs. 

2-5  to  be  a  measure  of  relative  error,  it  is  seen  to  be  on  the  order  of 
5  percent. 


C.  Discussion  and  Comparisons 

Comparisons  hove  been  made  between  the  empirical  curve  of 
Fig.  3,  found  for  the  experimental  values  of  Cs^-37,  concrete. 
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0089  -  0.50,  and  the  values  of  Clifford  and  of  Chilton -Huddleston. 

In  Fig.  11  the  data  of  Clifford  are  represented  by  point  symbols 
corresponding  to  particular  values  of  0.  It  may  be  seen  that,  though 
the  scatter  of  data  points  is  large,  the  data  are  not  inconsistent  with 
the  empirical  curve.  The  values  of  Chilton -Huddleston  were  converted 
to  and  are  shown  in  Tig.  12  together  with  the  empirical  curve.  Here 
values  derived  from  the  C-H  formula  vary  smoothly  with  9s  but  there 
appear  to  be  systematic  variations  with  both  9  and  9g  from  the  experi¬ 
mentally  based  curve . 

The  latter  conclusions  were  verified  when  a  point  by  point 
comparison  was  made  by  calculating  the  ratio  of  the  NRDL  data,  presented 
in  Tables  1  and  2,  to  the  C-H  formula  values  using  the  Raso  parameters, 
presented  in  Tables  3  and  It  was  found  that  agreement  is  less  than 
satisfactory.  The  difference  between  theory  and  experiment  is  rarely 
less  than  20  percent  as  shown  by  the  histogram  for  Cal37  in  Fig.  13  and 
frequently  exceeds  50  percent  as  shown  by  the  histogram  for  Co®°  in  the 
same  figure. 

The  original  Monte  Carlo  calculations  of  Raso  for  concrete  dose 
albedo  have  been  redone  following  discovery  of  an  error  In  the  program. 
Chilton-Huddleston  have  derived  new  values  for  the  parameters  In  their 
semi -empirical  formula  based  on  the  revised  data  of  Raso.  The  new 
values  for  08*37  are:  C  -  0.04o4  ±  0.0020,  C*  -  0.0172  i  0.0012;  and 
for  Co®0  are  €  -  0.0645  i  0.0022,  0*  »  0.0090  1  0.0007*  No  tabulation 
of  the  dh  values,  based  on  these  new  parametric  values,  was  available 
at  the  time  this  report  was  compiled.  Therefore,  the  %  values,  based 
on  the  original  parametric  values  noted  in  Tables  3  and  4,  have  been 
tabulated.  The  large  discrepancy  noted  above  in  the  point  by  point 
comparisons  cannot  be  explained  by  these  parametric  differences  however, 
since  they  would  increase  the  new  theoretical  values  of  <%  by  no  more 
than  3  or  4  percent. 

Although  it  is  possible  that  the  discrepancies  found  between 
these  results  and  the  Monte  Carlo  calculations  are  due  to  some 
undetected  systematic  error,  we  feel  it  is  more  likely  that  they  ax's 
due  to  the  use  of  the  Klein-Nlshina  cross  sections.  Considering  only 
the  incoherent  scattering  from  free  electrons  is  an  oversimplification 
which  is  not  physically  accurate.  While  most  of  the  neglected  processes 
contribute  only  a  small  amount  to  penetration  data  it  is  expected® 
that  these  effects  might  show  up  in  comparisons  between  theoretical  and 
experimental  albedo  values. 


*  Private  eoammleatlon. 
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Fig.  11  Comparison  of  Clifford  with  NRDL;  A^(ft)  vs.  08.  The  experi¬ 
mental  data  of  Clifford®  are  plotted  with  symbols  corresponding  to  the 
0  angles  at  which  they  were  taken. 
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8,  (RADIANS) 


Fig.  12  Comparison  of  Chilton-Huddleston  with  KRDLj  A^(fi)  vs.  es.  ! 
calculated  data  of  Chilton-Huddleston  are  plotted  with  symbols  corre 
sponding  to  the  6  angles  for  which  they  were  calculated. 


FREQUENCY 


ad(&)  NRDL 

RATIO:  - - - 

Qd  CHILTON -HUDDLESTON  FORMULA 


Fig.  13  Histograms  of  the  Comparison  Ratio  Distributions, 
ad(n)  (NRDL) 

ad  (Chilton -Huddleston  formula) 
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APPENDIX  I 


DISCUSSION  07  AIBEDO 


A.  General 


Literally,  albedo  means  whiteness;  technically,  it  Is  a 
reflection  coefficient.  It  has  been  appropriated,  from  its  original 
usage  in  astronomy  and.  optics  to  help  describe  gamma-ray  scattering. 
Since,  in  contrast  to  light  photons,  gamma  photons  interact  throughout 
a  volume  rather  than  at  a  surface,  the  concept  of  gamma  albedo  needs  to 
be  examined  closely. 


B.  Notation 

o  Jj. 

We  follow  the  notations  of  Goldstein0  and  Raso  where  possible. 
The  current-ref lection  coefficient  (albedo)  is  always  represented  by  the 
symbol  a,  the  flux-reflection  coefficient  by  A,  flux  density  by  I  and 
current  density  by  J.  The  subscripts  n  and  d  denote  number  and  dose, 
respectively.  These  and  other  symbols  used  in  this  report  are  listed 
alphabetically  in  the  glossary. 

Degrees  of  differentiation  of  a  functional  quantity  will  be 
represented  by  the  number  of  variables  appearing  as  arguments  of  the 
function.  Thus:  C^(fl,ho)  in  differential  angular  energy  dose  albedo. 

C.  Definitions 

The  fundamental  physical  quantity  in  gamna-ray  transport  is 
the  flow  of  photons.  Specifically: 


1-1  In(r,iJ,ho)  I  photons  «cm”2.s”1*8ter-1-kev~1 


is  the  differential  angular  energy  distribution  of  the  photon  mober 
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flux  density  and,  when  multiplied  by  dftl(ho),  will  give  the  number  of 
photons/s  crossing  a  unit  area  at  r  whose  normal  lies  along  the  unit 
direction  vector  n,  the  phgtons*  energies  lying  between  ho  and  hvH-d(ho) 
and  their  directions  along  fl  within  the  solid  angle  dfl.  See  Pig.  14. 

For  the  description  of  albedos  or  reflection  coefficients  it  is 
convenient  to  set  up  the  coordinate  axes  so  that  the  reflecting  surface 
is  the  xy  plane  and  so  that  the  point  of  intersection  of  the  source 
collimator  axis  with  this  plane  determines  the  origin. 

For  this  experiment  the  experimental  quantities  are  independent 
of  x  and  y  and  z  s  0  since  we  have  placed  the  origin^of  our  coordinate 
system  at  the  surface.  We  thus  drop  r  and  write  3^(0, ho) . 

If  the  reference  unit  area  at  the  origin  remains  in  the  back* 
scattering  surface  as  fl  varies,  we  have  a  different  quantity  called 
angular  energy  number  current  density,  which  is  related  to  angular 
energy  flux  density  by: 


1-2  Jn(n,ho)  »  In(iJ,ho)  cose 


where  9  is  the  angle  between  0  and  the  surface  normal. 


D.  Integral  Albedo 


JL  o) 


If,  in  a  time  t,  N'  incident  gamma  photons  penetrate  a 
surface  area,  a,  of  a  medium;  and  if,  after  scattering  in  the  medium, 
N  photons  come  back  out  through  the  same  area;  then,  the  Integral 
gamma-ray  number  albedo  for  these  incident  photons  and  this  medium  is 
defined  to  be: 


1-3 


°n  *  »/* 


(°) 


Note  that  in  this  example  the  incident  current  density, 
i  Nv°)  a-lt-1  and  the  reflected  current  density  is  Just  Jn  « 

N  a_lfi.  Thus  we  see  that  albedo  may  be  defined  as  a  ratio  of  current 
densities  since: 


Jn(o) 
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1-4 


Jn  N/at 


N _  „ 


E.  Differential  Albedo 


One  could  define  many  differential  albedos,  each  one  corres¬ 
ponding  to  different  combinations  of  incident  and  scattered  currents  of 
number,  energy  or  dose,  each  being  partially  or  totally  differentiated 
with  respect  to  direction  and  energy.  Of  these  many  possibilities  the 
albedo  quantity  reported  in  this  experiment  is  the  angular  dose  albedo 
for  an  "infinitely"  thick*  slab  irradiated  by  an  "infinitely"  broad** 
beam  of  approximately  parallel  gamma  rays.  Symbolically: 

oo 


/ 


d(ho) 


1-5 


ad(a) 


n  n°° 

J  J  J  Jd^O\eS0,,hb,)sin»*d©*d^*d(hu) 


where: 


Jd(M,ko)  -  Id(e,0,hu)coB®;  ld(e,0,hv>)  *  ho 


na(hv>) 

-2- - In(o,0,ht>). 


is  the  air  energy  mass -absorption  coefficient.) 


*  Bo  effect  is  observed  when  additional  slab  material  la  added. 

**  Bo  effect  is  observed  when  the  beam  radius  is  increased. 

***  Id(D)  =  /  ld(©,0,hu)  d(hu)  is  the  differential  angular  first 

o 

collision  air  absorbed  "dose”  rate;  referred  to  as  angular  dose 
rate. 
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Since,  in  this  experiment,  the  source  gamma  rays  were  both  mono- 
directional  and  monoenergetic. 


Jd(o)(e',0«,ho‘)  =  (Jd)o  B(9'-e0)  6(0«)  6(hp»-hoo) 


where 

(Jd>0  ■  Ja<o)(eo'°'huo>- 


Integration  of  1-5  then  yields 


1-6  ad(n)  =  jjj-  . 

'  d'o 

4 

This  differential  albedo  is  identical  to  Raso’s  a. 
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APPENDIX  II 


EXPERIMENTAL  APPARATUS  AND  METHODS 


A.  Experimental  Design 

The  quantity  to  be  determined,  o^(0),  imposes  several 
conditions  upon  the  experimental  apparatus  and  method.  First,  the  solid 
angle  (at  the  center  of  the  slab  face)  subtended  by  the  dosimeter  should 
be  small  since  the  experimental  results  are  to  be  differential  in 
scattering  angle.  Second,  the  slab  area  "seen"  by  the  dosimeter  should 
be  well  defined  since  this  area  is  used  to  convert  dose  rate  to  "dose- 
current."  Uiird,  the  slab  area  and  thickness  should  be  "infinitely" 
large,  as  noted  in  Appendix  I.  Fourth,  the  slab  should  be  uniformly 
irradiated.  The  experimental  configuration  chosen  to  satisfy  these 
conditions  is  shown  in  Fig.  15a. 

Specifically,  the  first  condition  was  met  by  using  a  lead- 
shielded  plastic  scintillation  detector  with  a  collimator  in  the  form 
of  a  long  cylindrical  hole  in  the  detector  shield.  The  solid  angle 
subtended  at  the  center  of  the  slab  face  by  the  exposed  area  of  the 
scintillator  has  been  determined  to  be  223  p-steradians  (see  Fig.  15b). 
Because  of  the  small  cross  section  for  backscattering,  the  extremely 
small  solid  angle  subtended  by  the  detector  and  available  source-strengths 
(Co°0,  — 11.5  curies;  Csl37,  ~"8  curies),  the  level  of  dose  rate  measured 
in  this  experiment  was  from  0.2  to  2.0  mr/hr.  To  accomplish  this  measure¬ 
ment  with  some  degree  of  reliability  in  a  reasonable  length  of  time  it 
was  necessary  to  use  a  scintillation  dosimetry  system  similar  to  that 
used  by  Dahlstroin.3 

Although  the  second  condition  could  be  met  by  increasing  the 
(collimator  length )/( detector  distance)  ratio,  the  resulting  detector- 
shield  weight  and  interference  with  the  incident  beam  make  this  solution 
impractical.  The  relationships  between  the  irradiated  area  of  the  slab 
and  the  effective  "viewed"  area  are  shown  in  Fig.  15c  and  are  discussed 
in  detail  in  reference  3- 

While  the  umbral  area  is  only  0.44  secein  the  penumbral  2 
region  increases  the  effective  "viewed"  area  (a  sec©)  to  10.78  sec©in^ 
and  thus  the  detection  solid  angle  (H  =a  /R2)  to  5330  p-steradians . 
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APPENDIX  II 


EXPERIMENTAL  APPARATUS  AND  METHODS 


A.  Experimental  Design 

The  quantity  to  be  determined,  Q^(0),  imposes  several 
conditions  upon  the  experimental  apparatus  and  method.  First,  the  solid 
angle  (at  the  center  of  the  slab  face)  subtended  by  the  dosimeter  should 
be  small  since  the  experimental  results  are  to  be  differential  in 
scattering  angle.  Second,  the  slab  area  "seen"  by  the  dosimeter  should 
be  well  defined  since  this  area  is  used  to  convert  dose  rate  to  "dose- 
current."  Uiird,  the  slab  area  and  thickness  should  be  "infinitely" 
large,  as  noted  in  Appendix  I.  Fourth,  the  slab  should  be  uniformly 
irradiated.  The  experimental  configuration  chosen  to  satisfy  these 
conditions  is  shown  in  Fig.  15a. 

Specifically,  the  first  condition  was  met  by  using  a  lead- 
shielded  plastic  scintillation  detector  with  a  collimator  in  the  form 
of  a  long  cylindrical  hole  in  the  detector  shield.  The  solid  angle 
subtended  at  the  center  of  the  slab  face  by  the  exposed  area  of  the 
scintillator  has  been  determined  to  be  223  p-steradians  (see  Fig.  15b). 
Because  of  the  small  cross  section  for  backscattering,  the  extremely 
small  solid  angle  subtended  by  the  detector  and  available  source-strengths 
(Co°°,  ''■'11.5  curies;  Csl37,  ^8  curies),  the  level  of  dose  rate  measured 
in  this  experiment  was  from  0.2  to  2.0  mr/hr.  To  accomplish  this  measure¬ 
ment  with  some  degree  of  reliability  in  a  reasonable  length  of  time  it 
was  necessary  to  use  a  scintillation  dosimetry  system  similar  to  that 
used  by  Dahlstrom.3 

Although  the  second  condition  could  be  met  by  Increasing  the 
(collimator  length )/( detector  distance)  ratio,  the  resulting  detector- 
shield  weight  and  interference  with  the  incident  beam  make  this  solution 
impractical.  The  relationships  between  the  irradiated  area  of  the  slab 
and  the  effective  "viewed"  area  are  shown  in  Fig.  15c  and  are  discussed 
in  detail  in  reference  3. 

2 

While  the  umbral  area  is  only  0.44  sec©in  the  penumbral  2 
region  increases  the  effective  "viewed"  area  (a  sec©)  to  10.78  sec©in 
and  thus  the  detection  solid  angle  (fl  =a  /R2)  to  5330  p-steradians. 
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1T7 

Fig.  15a  Experimental  Configuration.  The  Cs  source  is  positioned 
for  an  incidence  angle  of  arcos  0.50.  The  detector  is  positioned  for 
6  =  45°  and  0  =  112.5°.  3he  slab  material  is  concrete.  The  Co^O 
source  is  shown  stored  and  secured. 
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The  third  condition,  implied  by  specifying  the  slab  geometry 
as  semi -infinite,  means  that  any  increase  in  the  irradiated  slab  area 
or  the  slab  thickness  must  not  result  in  a  change  in  the  albedo  from  the 
viewed  area.  A  minimum  irradiated  circular  area  of  15  in.  radius  and 
slab  thicknesses  of  4  mean  free  paths  or  greater,  coupled  with  a  maximum 
"viewed”  area  having  a  semi -major  axis  of  6  in.  were  chosen  to  satisfy 
this  condition. 9 

By  "uniformly  irradiated"  in  condition  four  we  mean  that  the 
incident  gamma  rays  are  essentially  monoenergetlc,  monodirectional  and 
of  equal  flux  density  over  the  irradiated  area.  Condition  four  was 
approximately  satisfied  by  placing  collimated  Cs^-37  and  Co^°  "point" 
isotropic  sources  at  a  distance  of  97  in.  from  the  slab.°  The  dose 
rates,  measured  by  calibrated  Landsverk  and  Baldwin-Farmer  dosimeters, 
under  these  conditions  were  5.25  x  10^  and  2.11  x  lo7  keV/gm  sec  for 
Cs1^  and  Co^O  respectively.  These  quantities,  adjusted  for  decay 
during  the  course  of  the  experiment,  are  the  incident  dose  rates,  (i.)  , 
used  in  calculating  a.(  it) .  0 


B.  Dosimetry  System 

The  essential  details  concerning  the  detection  system  used  in 
this  experiment,  the  calibration  of  the  system,  and  the  means  by  which 
the  output  of  the  system  is  converted  into  detected  dose  rate  is 
reported  here.  First,  the  sections  shown  in  Fig.  l6  are  discussed, 
followed  by  preoperation  calibration,  dosimetry  system  calibration,  and 
data  acquisition  and  recording. 

1.  Detector. 

The  gamma  sensor  is  a  1  in.  diameter  by  a  1  in.  height 
plastic  scintillator  of  a  polymerized  solution  of  p-terphenyl  and 
tetraphenylbutadiene  in  styrene. 

The  thickness  of  the  lead  shielding  around  the  plastic 
scintillator  varied  from  10  in.  in  the  forward  direction  to  3  in.  at  the 
sides.  A  3/4  in.  diameter  by  10  in.  long  aperture  in  the  shielding 
acted  as  a  collimator  to  define  both  the  exposed  frontal  area  of  the 
scintillator  and  the  "viewing"  solid  angle  of  the  detector. 

The  light  from  this  scintillator  is  converted  into  current 
pulses  by  a  l4-stage  SMI-95 3&S  photomultiplier  operated  at  1100  volts. 
The  photomultiplier  output  is  coupled  to  a  two-stage  Franklin  pre¬ 
amplifier.  Hie  first  stage,  a  one  tube  class  A  amplifier,  is  DC 
coupled  to  a  cathode  follower.  The  overall  preamplifier  gain  is  3  with 
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Fig.  1 6  Dosimetry  System  Diagram. 


sufficient  negative  feedback  to  provide  excellent  gain  stability. 

The  output  impedance  of  the  preamplifier  is  low  to  match 
the  impedance  of  the  25  ft  length  of  coaxial  cable  connecting  the 
detector  and  the  converter  sections . 

2.  Converter. 

The  converter  section  is  composed  of  a  Hamner  Model  N  303 
non-overloading,  linear  amplifier,  the  analog  to  digital  converter  unit 
(ADC)  from  a  Penco  Model  PA-4  pulse-height  analyzer  and  a  NRDL  Model  1 
Calibration  and  Control  Unit.  An  upper  bound  discriminator  was  added 
to  the  ADC  to  reject  pulses  exceeding  the  upper  bound  voltage,  V  . 

The  amplifier  gain  was  standardized  to  spread  the  pulse- 
height  distribution  from  350  keV  photons  over  the  full  100  volt  ADC 
range.  (See  Fig.  17*) 

With  this  standard  system  amplification  some  pulses 
exceed  the  upper  bound  voltage  for  photon  energies  above  350  keV.  By 
adjusting  the  amplifier  gain  to  one  half  its  standard  value  and  raising 
the  lower  ADC  discriminator  level  to  l/2  Vu  it  is  possible  to  convert 
only  these  excess  pulses  for  photon  energies  up  to  700  keV.  By  this 
two  gain  technique  the  dosimetry  system  is  able  to  handle  the  range  of 
photon  energies  of  interest  in  this  experiment. 

Hie  ADC  unit  will  convert  any  analog  pulse  greater  than 
the  minimum  threshold,  =  4  volts  and  less  than  the  upper  bound 
voltage,  Vu  =  104.5  volts.  For  each  conversion  process  an  events  pulse 
is  generated.  In  the  amplitude  to  time  conversion  process  a  one 
megacycle/s  pulse  generator  is  gated  on  for  a  length  of  time  related  to 
the  peak  voltage  amplitude,  V,  attained  by  the  analog  pulse.  Hie 
relationship  is  such  that  the  number  of  digital  pulses,-  q,  in  this  pulse 
train  is  numerically  equal  to  the  quantity  (V  -  Vrp)  where  the  bar 
indicates  that  V  -  Vrp  is  rounded  off  to  the  nearest  integer  voltage. 

Thus,  the  presence  of  the  threshold  produces  aberrations  in 
the  response  in  two  ways :  No  pulse  less  than  4  volts  in  height  is 
converted;  and  each  q  is  4  units  smaller  than  the  corresponding  analog 
pulse  height.  Corrections  for  these  threshold  effects  have  been 
incorporated  into  the  dose -rate  formula  of  the  reduct or  and  are 
discussed  in  Section  D  (Dosimetry  System  Calibration). 

3 .  Register . 

The  digital  register  accumulates  the  digital  pulses,  q, 
generated  by  the  converter.  The  total  number  of  such  pulses  accumulated 


44 


PULSE  HEIGHT  (keV  UNITS) 


Fig.  17  Pulse -Height  Distribution.  Ihe  distribution  is  of  the  response 
of  the  1  in.  diameter  by  1  in.  height  plastic  scintillator,  used  in  the 
experiment,  to  Csl37  gamma  rays . 


*5 


in  a  time  t  is  designated  by  Q  =  £q. 


The  events  register  accumulates  the  events,  P,  converted  in 

time  t. 

The  timing  register  accumulates  the  total  number  of  kc/s 
pulses  generated  by  the  live-timing  circuit.  Ihus,  the  timing  register 
divided  by  a  thousand  gives  the  elapsed  counting  time,  t,  in  seconds. 

4.  Reductor. 

In  the  reductor  the  registered  quantities  Q,  p  and  t  are 
combined  to  give  detected  dose  rate  by  the  formula: 


'  total 


■) 


^total  ~  %  +  2%/2  8111  Ptotal  Ps  +  2Ps/2 


where:  Qg  and  Ps  are,  respectively,  the  net  registered  digital  pulses 
and  the  net  registered  events  for  the  standard  mode;  and  Qg/g  and  Ps/2 
are,  respectively,  the  net  registered  digital  pulses  and  the  net  ' 
registered  events  above  the  standard  mode  upper  bound  taken  with  half¬ 
standard  gain.  Kg  is  an  empirically  determined  dose-rate  conversion 
factor  and  T  is  a  threshold  correction  parameter. 

It  should  be  noted  that  the  values  of  Kg  and  T  are  of 
critical  importance  in  this  experiment  since  they  convert  the  system 
into  a  dosimeter.  The  determination  of  these  values  will  therefore  be 
discussed  at  some  length.  However,  prior  to  this  discussion,  it  is 
desirable  to  dispense  with  the  details  associated  with  the  preoperation 
calibration  of  the  electronics . 

In  addition  to  the  above  sections  a  calibration,  control 
and  timing  unit  serves  to  make  the  dosimetry  system  operable.  It 
provides  a  method  for  calibrating  the  system,  a  preset  live-time 
capability  and  certain  control  functions  not  found  in  the  PA-4  ADC  unit. 


C.  Preoperfction  Calibration 
1 .  ADC  Calibration . 
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Reference  voltage,  channel  linearity,  channel  width  and 
threshold  checks  of  the  ADC  unit  vere  made  as  specified  In  the  Penco 
PA-4  Instruction  Manual. 

She  upper  hounds  circuit  of  the  ADC  vas  adjusted  to  prevent 
the  conversion  of  any  analog  pulse  -which  would  give  a  digital  pulse 
train  value,  q,  greater  than  100.  Since  the  range  of  analog  pulse 
voltage  amplitudes  which  are  converted  to  a  q  of  100  Is  from  103.?  v  to 
104.5  v,  this  adjustment  was  performed  while  converting  104.5  v  pulses 
from  external  60  c/s  pulse  generator.  Because  there  is  a  symmetrical 
distribution  of  pulses  around  the  average  value  of  104.5  v,  half  of  these 
pulses  should  he  converted  if  the  upper  hound  discriminator  Is  working 
properly.  Thus  an  events  scaler  rate  p/t  of  30  s“l  vas  the  objective  In 
adjusting  the  upper  bound  discriminator. 

2.  Linearity  and  Threshold. 

Upon  completion  of  the  ADC  calibration,  an  Investigation  of 
the  linearity  and  threshold  of  the  system  from  the  preamplifier  through 
the  ADC  vas  made.  Ten  second  sampling  periods  vere  used  to  measure  Q 
and  P  values  for  various  pulser  settings .  A  plot  of  q  ■  <^P  versus 
pulser  setting  resulted  in  a  linearity  curve  with  the  threshold  being 
the  Intercept.  Figure  l8  is  a  plot  of  typical  results. 


D.  Dosimetry  System  Calibration 

Ideally,  the  accumulated  digital  counts,  Qtotal>  divided  by  the 
counting  time  t  is  proportional  to  dll,  However,  as  mentioned  in  the 
converter  discussion,  the  threshold  of  the  ADC  unit  introduces  two  types 
of  aberration.  The  threshold  subtractive  effect  reduces  each  analog 
pulse  voltage  by  4  volts  and  the  threshold  loss  effect  Ignores  a 
fraction  of  the  input  pulses .  This  fraction  Increases  as  the  photon 
energy,  ho,  decreases.  In  addition,  the  fraction  of  energy  absorbed  in 
the  plastic  scintillator  due  to  multiple  collisions,  rather  than  first 
collisions,  also  Increases  as  the  photon  energy,  ho,  decreases.  The 
formula, 


dD  - 


of  the  reductor  section  is  an  attempt  to  compensate  for  these  energy 
dependent  aberrations.  The  method  of  determining  K  by  adjusting  T,  to 
minimize  these  aberrations,  follows. 


q  (DIGITAL  PULSES/ EVENTS  PULSE) 


NRDL  379-64 


Fig.  l8  Linearity  Calibration  Curve. 
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Four  calibration  sources  were  used,  each  producing  monoenergetic 
gamma  rays.  Each  source  was  placed  on  the  collimator  axis  as  shown  in 
Fig.  19.  The  dosimeter  measurements 


were  taken  for  each  calibration  gamma  photon  energy,  hi^. 

Without  disturbing  the  source  and  detector  shield  .geometry  the 
plastic  scintillator  and  photomultiplier  were  replaced  by  a  1  in.  by 
1  in.  diameter  Nal(Tl)  scintillator  and  photomultiplier.  The  amplified 
output  pulses  from  this  detector  were  analyzed  and  recorded  by  a  Penco 
PA-3  analyzer.  From  these  pulse -height  distributions  and  the  known 
detection  efficiency  and  detection  area  of  the  Nal  detector,  the  incident 
photon-number-flux  density,  In(ho)d,  was  calculated  for  each  source. 

Hie  dose  rate  at  the  detector  was  calculated  by  the  formula: 


(ho)  In(ho) 


i‘ 


Dose -rate  conversion  factors  for  each  calibration  gamma  photon 
energy,  ho^,  were  calculated  by  the  formula: 


K(ho,T)t  = 


Id(hb)/ 


+  T 


total 


i 


for  several  values  of  T. 

Figure  20  shows  the  energy  dependence  of  K(ho, T)  for  four 
values  of  T.  T  =  10  was  the  final  value  chosen  to  calculate  each 
K(ho)i.  The  system  dose-rate  conversion  factor.  Kg,  was  obtained  by 
taking  a  simple  average  of  the  four  K(hp,  10^  factors: 

Ks  -  l/4  K(hb,10)i  =  0.440  ±  0.006. 

i=l 
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Fig.  19  Dosimetry  System  Calibration  Geometry. 
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Fig.  20  Energy  Dependence  of  Dose  Conversion  Factor.  The  set  of  points 
marked  T  =  0  shows  the  uncorrected  factor.  The  points  labeled  T  =  4  are 
corrected  for  the  threshold  subtractive  aberration.  The  points  for  T  = 
10,  the  value  finally  adopted  compensates  for  the  threshold  loss 
aberration.  The  points  for  T  =  12  shows  the  effect  of  an  additional 
compensation.  The  solid  line  labeled  Kg  is  an  average  of  the  four  T  * 
10  points  and  is  the  system  dose  conversion  factor  used  in  formula  of 
the  reductor  section. 
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Inserting  these  values  of  Kg  and  T  into  the  formula  of  the 
reductor  section  we  have: 


dD  =  o.44o 


^  Ssotal 


+  10 


‘  total ^ 


keV  gm""1'  s” 


l]- 


E.  Data  Acquisition  and  Recording 

The  variable  parameters  were  the  radiation  sources,  Cs1^  and 
Co°  ;  the  slab  materials,  concrete,  aluminum  and  steel;  the  incidence 
angles,  0O  =  cos  1.00,  0.75,  0.50;  and  the  detection  angles,  6  =?  8°, 

15°,  22°,  32°,  1+5°,  60°  and  72°  and  0  =  7.5°,  22.5°,  37-5°,  52.5°,  67.5°, 
82.5°,  97- 5°,  112.5°,  127. 5°,  142.5°,  157-5°  and  172. 5°. 

An  almost  complete  set  of  data  for  all  detection  angles  was 
taken  for  the  concrete,  but  for  the  aluminum  and  steel,  measurements 
were  taken  for  selected  values  of  9  and  0. 

The  slabs  were  exposed  and  data  were  taken  for  sampling  times 
of  1000  seconds.  For  each  combination  of  the  variable  parameters 
described  above  a  foreground  and  background  was  made  at  standard  gain 
and,  where  necessary,  at  half  standard  gain.  The  background  runs  were 
made  with  the  source  "on"  but  with  a  lead  plug  in  the  detector  collimator. 

The  net  readings  (Sforeground-Sjackground)  and  (^foreground- 
Dackground)  were  combined  to  give  the  dose  rate  at  the  detector. 

Daily  checks  were  made  of  the  system  amplification  by  feeding 
calibration  pulses  of  a  standard  amplitude  into  the  preamplifier  in  the 
same  way  as  for  the  linearity  and  threshold  calibration.  Any  shift 
noted  was  compensated  for  by  adjusting  the  amplifier  fine  gain  control. 
Following_this  adjustment,  the  overall  system  gain  was  measured  by 
making  a  q.  measurement  for  a  Hg20-^  calibration  source.  The  distribution 
of  the  daily  Hg2®3  q  values  had  a  maximum  variation  from  the  average  of 
less  than  1  percent  during  the  period  of  data  taking. 
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